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Abstract 

We present a study of the one dimensional Su-Schrieffer-Heeger model Hamil- 

tonian by a diagrammatic perturbative method in the weak electron-phonon 

coupling regime. Exact computation of both the charge carrier effective mass 

and the electron spectral function shows that electrons are good quasiparti- 

cles in the adiabatic and antiadiabatic limits but novel features emerge in the 

intermediate regime, where the phonons and the electrons compare on the 

energy scale. Together with a sizeable mass enhancement we observe, in the 

latter regime, a spread of the spectral weight (among several transition peaks) 

associated with an increased relevance of multiphonons contributions at larger 

e-ph couplings. Accordingly electrons cease to be the good quasiparticles and 

an onset of polar on formation is favoured. 
PACS: 63.20.Dj, 71.18.+y, 71.38.+i 



1 



Typeset using REVTgX 



l.INTRODUCTION 



The coupling of the electrons to the quantum lattice vibrations is a key issue in solid 
state physics. While in the weak coupling limit electrons generally behave as good quasipar- 
ticles, in strong coupling regimes electrons may drag the lattice deformation induced by the 
scattering process. In the latter case, the charge carrier is the unit comprising electron and 
dressing cloud of phonons, namely the polaron [ED]. A wide literature has been produced on 
the subject during the last decades [0-171] and particular emphasis has been recently laid on 
the polaronic properties of high T c superconductors ||. Theoretical investigations usually 
start from the Holstein molecular crystal model f| which assumes a momentum independent 
coupling of electrons to dispersive optical phonons leading to polaron formation if the en- 
ergy gain associated with the lattice deformation is larger than the kinetic energy due to the 



tight binding hopping in the bare band |10|. As a remarkable feature of the polaron physics 
|rTHl4[ , the strength of the coupling drives the transition from a large size polaron (at weak 
couplings) to a small size polaron (at strong couplings). This transition []15j is accompanied 
by a sizeable enhancement of the effective mass whose value may change considerably ac- 
cording to the degree of adiabaticity and the peculiarities of the lattice structure |16| , |17j . In 
some systems however |jl8|-p0f the e-ph interaction modifies the electron hopping matrix ele- 
ments thus leading to a momentum dependent coupling function. In this case an appropriate 



theoretical framework is offered by the Su-Schrieffer-Heeger(SSH) model Hamiltonian |2T 



originally proposed to explain the conducting properties of quasi one dimensional polymers 
as polyacetylene f22| whose twofold degenerate ground state sustains well known nonlinear 
local excitations, the solitons. More generally, the SSH Hamiltonian provides an alternative 
(to the Holstein model) tool to analyse the physics of polaron formation as tuned by the 
strength of the e-ph coupling in one dimensional systems. This physics has been less deeply 
investigated than in the Holstein model. With the present paper we address precisely this 
question attacking the SSH Hamiltonian by a weak coupling perturbative method. Although 
this approach is clearly not adequate to capture the multiphononic nature of the polaronic 
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quasiparticle, still it can provide useful informations regarding the onset of polaron forma- 
tion in some portions of parameter space. Here we look first at the mass renormalization 
exploring a wide range of values for the adiabatic parameter and, successively, we compute 
the electronic spectral function to detect whether and to which extent bare electrons behave 
as good quasiparticles. The Section 2 outlines the SSH model and contains the results of 
this study while some conclusions are drawn in Section 3. 



2. MODEL AND RESULTS 



In ID, the real space SSH Hamiltonian reads 



H= E Jr,r + l(frfr+l + fl+lfr) + E(|^ + f K " ^r+lf) 

Jr,r+1 = ~[J + «(«r - Ur+l)] (0-1) 

where J is the nearest neighbors hopping integral for an undistorted chain, a is the 
electron-phonon coupling, u r is the dimerization coordinate which specifies the displacement 
of the r— lattice site from the equilibrium position along the molecular axis, p r is the 
momentum operator conjugate to « r , M is the ion (ionic group) mass, K is the effective 
spring constant, f} and f r create and destroy electrons on the r— site. Let's expand the 
lattice displacement and its conjugate momentum in terms of the phonon creation and 
annihilation operators b q and b q and Fourier transform the electron operators 

u r = E Tz n \ AT ( b -i + K) exp(igr) 
q J2MNu q 



fr= -~^T,*Mikr)f k (0.2) 



in order to obtain the SSH Hamiltonian in momentum space: 
H = H + H int 
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H =J2 £ kflfk + U 1 h \ h 1 
k q 

Hint = $>(* + q, k){bl q + b q )f k+q f k 

k,q 

e k = —Jcos(k) 



2 A K ■ 2/9< 

J n = 4— sm - 

q M y 2' 



%0i 

g{k + q,k) = . (sin(fc + q) - sin(fc)) (0.3) 

where N is the total number of lattice sites. The phonon dispersion relation is defined in 
the range q G [0, tt]. Assuming a reduced Brillouin zone (|g| < vr/2) the spectrum displays 
both an acoustic and an optical branch ]23| . The model contains three free parameters: the 



hopping integral J, the zone boundary frequency u n = 2JK/M which coincides with the 
zone center optical frequency in the reduced zone scheme, the coupling constant a 2 /AK. 
The full electron propagator in the Matsubara Green's functions formalism is defined as: 

G(k,r) = - E("l) n / dTi...dT n (T T f h {T)H int (T X ) ■ ■H int (r n )fl(0)) n (0.4) 

n=0 J ° X /0 

where (3 is the inverse temperature, T T is the time ordering operator, < ... >o indicates 
that thermodynamic averages are taken with respect to the unperturbed Hamiltonian and 
only different connected diagrams contribute to any order n. I have computed exactly 
the self-energy terms due to one phonon (n = 2 in eq.(4)) and two phonons (n = 4 in 
eq.(4)) scattering processes which determine the renormalized electron mass m e ff through 
the relations: 

meff _ 1 - dReE k (e)/de\ k=0 . e= _j 
m 1 + dReL k {e)lde k \ k =Q; e =-j 

where, ReE k (e) = ReT$\e) + i?eSf a) (e) + ReE ( k h) (e) + ReT l < k 2c) {e) is the frequency 
dependent real part of the retarded self-energy. There are three contributions due to different 



connected two-phonons diagrams [21|. Their effect is however confined to the intermediate 
regime in which uj^ is comparable to the electronic energy J. Here the two phonons diagrams 
enhance the effective mass by ~ 15% with respect to the one phonon result. Instead, in 
the fully adiabatic and antiadiabatic regimes the two phonons contributions (evaluated at 



the band bottom) are negligible. Hereafter, the displayed results depend on the very one 
phonon self-energy Y^\ie m ) term (e m = (2m + 1)tt//3 with m integer number) whose finite 
temperatures analytic expression is given by: 



n B (uj q ) + n F (-e k - q ) n B (u q ) + n F (e k - q ) 



(0.6) 



^k—q ^q &k—q ^q 

ub and n F are the Bose and Fermi occupation factors respectively. We take a narrow 
band system setting J = O.leV. Figure 1 shows a sizeable mass enhancement in the inter- 
mediate adiabatic regime with a pronounced spike at cu^ ~ \[2J mainly due to scattering by 
\q\ = 7r/2-phonons. This feature is related to the ID electron and phonon dispersion rela- 
tions and it may be partly suppressed in higher dimensionality. However the onset of a mass 
renormalization starting at uo n ~ J/2 and, more evidently, at u) n ~ J (together with the in- 
creased relevance of multiphonons contributions) signals that polaron formation is expected 
in this regime while no mass enhancement is obtained in the adiabatic and antiadiabatic 
limits. Let's look now at the spectral function defined by A(k, e) = —2ImG re t(k, e) to get 



more insight [^] into the suggestions proposed by the effective mass behavior. In terms 
of the retarded self-energy, obtained by eq.(6) through analytic continuation e m — > e + i5, 
A(k, e) reads: 

A(k, e) = 2rc5(e - e k - i?eS fc ( e )) + (-)2ImZ k (e) (Q ?) 

(e - e k - -ReS fc (e)J + ImE k (e) 

The first addendum in eq.(7) contributes when JmSfc(e) = 0. The band bottom (k = 0) 
spectral function has been computed in a number of representative cases summing the self- 
energy term (eq.(6)) over 6000 g— points in the Brillouin zone and searching the zeros of the 
5-function argument (eq.(7)) and Im'E k=0 (e) over 10 5 points of the e— axis. The sum rule 

de 



/OO He 
-A(k = 0, e) = 1 
-oo Z7T 



has to be numerically fulfilled in principle. Here however we are approximating the total self- 
energy by the one-phonon term (eq.(6)) which linearly depends on the free parameter a 2 /AK. 
As the e-ph coupling grows multiphonons terms become more relevant and our approximation 



becomes less accurate. Accordingly, deviations from the sum rule are expected as a measure 
of the loss of spectral weight associated with higher order self-energy effects. In this regard, 
the sum rule numerical analysis permits to define the range of a 2 /4K values within which 
the one phonon approximation is reliable. For any choice of input parameters we are able 
to estimate the intrinsic error of our physical model. 

Figures 2(a) and 2(b) illustrate that electrons are good quasiparticles in the fully antia- 
diabatic regime: for both values of the e-ph coupling there is a well resolved peak due to 
the S— function contribution and located at e = — 98.8meV (Fig.2(a)) and e = — 90AmeV 
(Fig.2(b)), respectively. 

The sum rule is well satisfied in Figure 2(a) while a slight loss of spectral weight (8%) 
occurs in Figure 2(b) (where the dimensionless effective coupling is a 2 /(AKJ) = 0.1) , 
being ^ OQ deA{k = 0, e)/2n = 0.92. Infact, by increasing the e-ph coupling fourth order 
diagrams acquire relevance and should contribute to the fulfillement of the sum rule. Figures 
3 deal with the intermediate regime u n = J. Electrons are still good excitations in the 
extremely weak coupling case (Fig.3(a)) with one well defined transition at e = — 100.9meV 
but, at larger couplings (Fig. 3(b)), the sum rule is far from being satisfied and a 40% loss 
of spectral weight is observed together with a strong reduction of the main peak height. 
Here we see the failure of the electronic quasiparticle picture and the possible onset of a 
polaronic state. Let's come to the adiabatic regime represented in Figures 4. In Figure 
4(a) the sum rule is satisfied and the spectral weight spreads mainly in a few peaks around 
the highest one located at the energy e = — lOOmeV with a significant tail up to energy 
levels of order e ~ — SOmeV. At larger couplings (Fig.4(b)) the tail is appreciable up to 
e ~ — 50meV and the loss of spectral weight is ~ 20%. A well defined transition is however 
still present at e = —108.2meV. The temperature induced effects on the spectral function 
of the extremely weak coupling system are shown in Figures 5(a) for the T = 100K case and 
5(b) for the T = 300fT case, respectively. Comparison with the ground state case (Figure 
4(a)) emphasizes the progressive broadening of the main peak in a quasi continuum of states 
which, above room temperature, tend to have similar probabilities. 
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3. FINAL REMARKS 



Electrons are expected to be good quasiparticle and natural excitations of a system in 
which the electron-phonon interaction is weak. We have studied the Su-Schrieffer-Heeger 
Hamiltonian by a low order perturbative method assuming a coupling energy a 2 /4K always 
smaller than the bare hopping integral J which sets the scale of the electron energies. The ef- 
fective mass of the charge carriers varies with the degree of adiabaticity of the system: while 
no renormalization is observed in the adiabatic and antiadiabatic limits, a sizeable enhance- 
ment takes place when the phonons compare in energy with the electrons, namely in the 
intermediate regime. Although the pronounced mass peak may be partly reduced in higher 
dimensionality exact computations of fourth order diagrams clearly show that multiphonons 
effects become relevant in the intermediate regime thus signalling that a polaronic behavior 
may here occur. Analysis of the electronic spectral function corroborates these findings. 
More precisely than the effective mass, the electron spectral function, directly measures to 
which extent electrons are the good excitations of the system in any regime tuned by the 
adiabaticity parameter. Only one phonon scattering processes have been considered in the 
presented results. The spectral function sum rule has been assumed as a benchmark for the 
reliability of the one phonon approximation and the computed loss of spectral weight have 
been associated with an increased relevance of multiphonons contributions as a function of 
the strength of the e-ph coupling. 

While, in the antiadiabatic regime, the system presents well defined electronic transitions 
again, in the intermediate regime, electrons cease to be good quasiparticles as soon as the 
e-ph coupling grows, although well within a weak coupling scheme. This is an independent 
confirmation of the onset of polaronic features suggested by the mass calculation. Also in 
the moderately adiabatic regime we find that the spectral weight spreads over a region of 
energy levels but a defined transition is present in the ground state of the system. The 
polaronic properties of the model will be next studied in the strong coupling regime thus 
permitting a comparison with the Holstein molecular crystal model. Dimensionality effects 
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will be also investigated. 
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FIGURES 



FIG. 1. Renormalized masses (in units of bare band electron mass) versus the adiabaticity 
parameter, m^ff is due to the one phonon self-energy correction. 



FIG. 2. Electron spectral function in antiadiabatic regime and (a) extremely weak e-ph cou- 
pling; (b) moderately weak e-ph coupling. 
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FIG. 3. Electron spectral function in intermediate regime and (a) extremely weak e-ph coupling; 
(b) moderately weak e-ph coupling. 



FIG. 4. Electron spectral function in adiabatic regime and (a) extremely weak e-ph coupling; 
(b) moderately weak e-ph coupling. 
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FIG. 5. Finite temperatures electron spectral function in adiabatic regime and extremely weak 
e-ph coupling; (a) T = 100K, (b) T = 300K. 
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